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We investigate semi-inclusive photon-hadron production in the color glass condensate (CGC) 
framework at RHIC and the LHC energies in proton-proton (pp) and proton-nucleus (pA) collisions. 
We calculate the coincidence probability for azimuthal correlation of pairs of photon-hadron and 
show that the away-side correlations have a double-peak or a single-peak structure depending on 
trigger particle selection and kinematics. This novel feature is unique for semi-inclusive photon- 
hadron production compared to a similar measurement for double inclusive dihadron production 
in pA collisions. We obtain necessary conditions between kinematics variables for the appearance 
" of a double-peak or a single peak structure for the away-side photon-hadron correlations in pp 

and pA collisions at forward rapidities and show that this feature is mainly controlled by the 
ratio zt = Pr adron /py hoton . Decorrelation of away-side photon-hadron production by increasing the 
■ energy, rapidity and density, and appearance of double-peak structure can be understood by QCD 

saturation physics. We also provide predictions for the ratio of single inclusive prompt photon to 
hadron production, and two-dimensional nuclear modification factor for the semi-inclusive photon- 
hadron pair production at RHIC and the LHC at forward rapidities. 
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I. INTRODUCTION 



i-C ' It is generally believed that a system of partons (gluons) at high energy (or small Bjorken-x) forms a new state 
of matter where the gluon distribution saturates [l|. Such a system is endowed with a new dynamical momentum 
I scale, the so-called saturation scale at which non-linear gluons recombination effects become as important as the gluon 
radiation. The colorglass condensate (CGC) approach has been proposed to study the physics of gluon saturation 
at small-x region 0, yj]. The CGC formalism is an effective perturbative QCD theory in which one systematically re- 
sums quantum corrections which are enhanced by large logarithms of 1 /x and also incorporates non-linear high gluon 
£SJ ■ density effects. In the CGC approach, the main features of particle production at high energy remain universal and 
J> ' are controlled by the saturation scale. This picture has been successfully applied to many QCD processes from HERA 
to RHIC [H and the LHC [34Tl|. In this paper, we will employ the CGC formalism and show that the semi-inclusive 
photon-hadron (7 — h) production processes in pA collisions, i.e., p + A — > 7 + /i + A, offer more interesting insights 
. to the dynamics of gluon saturation. 

Photons radiated in hard collisions not via hadronic decays are usually called prompt photon. There are advantages 
0^ to studying prompt photon production as compared to hadron production. It is theoretically cleaner; one avoids the 
difficulties involved with description of hadronization and possible initial-state-final-state interference effects which 
may be present for hadron production and it can be therefore used as a baseline to interpret jet-quenching phenomenon 
in heavy ion collisions. A detailed studies of Ref. [l2[ showed that prompt photon production in pA collisions at RHIC 
and the LHC at forward rapidities is a sensitive probe of the small-x physics and QCD gluon saturation. On the 
same line, the semi-inclusive prompt photon-hadron production in pA collisions, has also advantages over a similar 
production of dihadron. In particular, in dihadron production, higher number of Wilson lines, and the Weizsacker- 
Williams and the dipole gluon distributions are involved [131 ] while in the photon-hadron production cross-section, 
only dipole gluon distribution appears [l2l . Il4j which is both experimentally and theoretically well-known, see for 
example Refs.@,[ISi3. 

Two particle correlations in high-energy collisions have played significant role to reveal QCD novel phenomena [22I — 
[24|. In particular, the photon-hadron jet correlations have been a very powerful probe of the in-medium parton energy 



loss in high-energy heavy- ion collisions [24|, |25[ . It was suggested in Ref. [12[ that the correlation of the back-to-back 
photon-hadron pair production in high-energy pp and pA collisions can be used to probe the gluon saturation at small- 
x region and to study the physics of cold nuclear matter in dense region. However, the correlation defined in Ref. [l2j 
may dep end crucially on the so-called underlying event and might be rather challenging to measure. Moreover, in 
Ref. [12| the correlation was studied in a very limited kinematics, see Sec. IV. In this paper, for the first time we 
study the coincident probability for photon-hadron correlation at RHIC and the LHC in both pp and pA collisions. 
Dihadron azimuthal angle correlation was already measured by the coincidence probability at RHIC [23(. We show 
that the away-side correlations for a pair of photon-hadron obtained via the coincident probability have a double or 
a single peak structure depending on kinematics and whether the trigger particle is selected to be a prompt photon 
or hadron. We obtain kinematics conditions for appearance of a double or a single peak structure for the away-side 
photon-hadron correlations which can be verified by the upcoming experiments at RHIC and the LHC. This novel 
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feature is unique for prompt-hadron production in contrast to dihadron production [l3j where the trigger particle 
can only be a hadron and it was already observed at RHIC [23[ that the way-side correlation has only single peak 
structure. The asymmetric nature of photon-hadron production, and the fact that in semi-inclusive photon-hadron 
production, QCD and electromagnetic interaction are inextricably intertwined, make the azimuthal correlation of the 
produced photon-hadron very intriguing. 

We will also provide quantitative predictions for the (two-dimensional) nuclear modification factor for the semi- 
inclusive photon-hadron production in pA collisions, and the ratio of single inclusive prompt-photon to hadron pro- 
duction in pp and pA collisions, at RHIC and the LHC at forward rapidities. 

This paper is organized as follows; In Sec. II, we first provide a concise description of theoretical framework by 
introducing the main formulas for the calculation of the cross sections of semi-inclusive photon-hadron (Sec. II- 
A), single inclusive prompt photon (II-B) and single inclusive hadron (ITC) production within the CGC approach. 
In Sec. II-D, we describe how to compute the main ingredient of our formalism, namely dipole-target forward 
scattering amplitude via the running-coupling Balitsky-Kovchegov evolution equation [l5j . In Sec. Ill we introduce 
the observables that we are interested to compute and our numerical setup. In Sec. IV, we present our detailed results 
and predictions. We summarize our main results in Sec. VI. 



II. THEORETICAL FRAMEWORK 



A. Semi-inclusive prompt photon-hadron production in pp and pA collisions 

The cross section for production of a prompt photon and a quark with 4- momenta p 1 and I respectively in scattering 
of a on-shell quark with 4-momentum k on a dense target either proton (p) or nucleus (A) at the leading twist 
approximation in the CGC formalism is given by [l4 |. 
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where y/~S is the nucleon-nucleon center of mass energy and the light-cone fraction x q is the ratio of the incoming 
quark to nucleon energies, namely x q = k~ / \J S/2. The pseudo-rapidities of outgoing prompt photon rj 7 and quark 
T]h are defined via p~ — ^|e'' 7 and l~ = -^e Vh . The subscript T stands for the transverse component. The vector &t 
denotes the impact-parameter of interaction. The angle between the final-state quark and prompt photon is denoted 
by A(j) and defined via cos(A</>) = l T P ^ T . Throughout this paper, we only consider light hadron production, therefore 

at high transverse momentum (ignoring hadron mass), the rapidity and pseudo-rapidity is the same. Note that due 
to the assumption of collinear fragmentation of a quark into a hadron, the angle A(f> is then the angle between the 
produced photon and hadron, assuming that the rapidity of the parent parton and the fragmented hadron is the 
same. In Eq. ([1]), Np(pT,x g ) is the imaginary part of (quark-antiquark) dipole-target forward scattering amplitude 
which satisfies the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation [2(| [27j 
and has all the multiple scattering and small- a; evolution effects encoded (see Sec. II-D). 

In order to relate the above partonic production cross-section to proton-target collisions, one needs to convolute the 
partonic cross-section in Eq. JT]) with the quark and antiquark distribution functions of a proton and the quark-hadron 
fragmentation function: 

■T — = / dx q f q (x q ,Q 2 ) D h/q (z f ,Q 2 ), (2) 

d 2 bx d 2 pT 1 d 2 pT drj^ drjh Jz"p in z f J d 2 bx d 2 px 1 d 2 lx drj^ drjh 

where p\ is the transverse momentum of the produced hadron, and f q (x q ,Q 2 ) is the parton (quark) distribution 
function (PDF) of the incoming proton which depends on the light-cone momentum fraction x q and the hard scale 
Q. A summation over the quark and antiquark flavors in the above expression should be understood. The function 
D} l / q (zf 1 Q) is the quark-hadron fragmentation function (FF) where Zf is the ratio of energies of the produced hadron 
and quark. 

The light-cone momentum fraction x qi x q ,x g in Eqs. (|1I2|) are related to the transverse momenta and rapidities of 
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the produced hadron and prompt photon via (see appendix in Ref. fl2j|). 
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B. Single inclusive prompt photon production in pp and pA collisions 

The prompt photon cross section in the CGC framework can be readily obtained from Eq. ([TJ by integrating over 
the momenta of the final state quark. After some algebra, the single inclusive prompt photon production can be 
divided into two contributions of fragmentation and direct photon [12| : 
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where p T is the transverse momentum of the produced prompt photon, and D 1 / q (z,Q 2 ) is the leading order quark- 
photon fragmentation function (28j . In order to relate the partonic cross-section given by Eq. Q to prompt photon 
production in pA collisions, we convolute Eq. (|4]) with quark and antiquark distribution functions of the projectile 
proton [20j | . 
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where a summation over different quarks (antiquarks) flavors is implicit. The light-cone fraction variables X q ^ X g i Z 111 
Eq. (|4|5|) are related to the transverse momentum of the produced prompt photon and its rapidity [l2| , 
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C. Single inclusive hadron production in pp and pA collisions 

The cross section for single inclusive hadron production at leading twist approximation, in asy mmetric collisions 
such as pA or forward rapidity pp collisions at high energy, in the CGC formalism is given by (30l. [3l|. 
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(7) 

where the variables r\h and p T are the pseudo-rapidity and transverse momentum of the produced hadron. The 
functions f q , N F (A) an d Dh/ q m the above are defined the same as in Eq. ([2]). The indices q and g denote quarks and 
gluon, with a summation over different flavors being implicit. The first two terms in the above expersion correspond 
to elastic contribution, namely an incoming parton scattering elastically with the CGC target [301 ] . This incoming 
parton with initial zero transverse momentum picks up transverse momentum of order saturation scale after multiply 
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scattering on the target. There is also inelastic contribution to the cross-section denoted by ^ meIastlc cor resp onding 
to a high transverse momentum parton radiated from the incoming parton in the projectile wave function [3ll |32|. In 
this case, the projectile parton interacts with target with small transfer momentum exchanges, but this is enough to 
decohere the pre-existing high-px parton from the hadron wave function and release it as an on-shell particle. The 
high-pT partons in the projectile wave function arise due to DGLAP splitting of partons. It was shown in Ref. [32] that 
at very forward rapidities the inelastic contributions are less important compared to elastic ones while it is significant 
at midrapidity at high-energy scatterings. 

The longitudinal momentum fractions x\ and X2 are defined as follows, 

iFW^i *i = — 5 x 2 = Xl e~ 2 ^. (8) 
v S z 

One should note that the light-cone fraction variables defined above for the inclusive hadron production are different 
from the corresponding light-cone variables for single inclusive prompt photon Eq. ([6]) and semi-inclusive photon- 
hadron Eq. ([3]) production. 



D. Small-x evolution equation and the dipole forward scattering amplitude 



The main ingredient in the cross-section of semi-inclusive photon- hadron production in Eq. (JTJ), single inclusive 
prompt photon production in Eq. Q and single inclusive hadron production in Eq. ([7]) is the universal dipole forward 
scattering amplitude which incorporates small-x dynamics and can be calculated via the first-principle non-linear 
JIMWLK equations [26|, [27|. In Eqs. (|1I4I7[) . the amplitude Np (Na) is the two-dimensional Fourier transformed of 
the imaginary part of the forward dipole-target scattering amplitude N"a(f) m the fundamental (F) or adjoint (A) 
representation, 



N A(F) (x,k T ) = / d 2 re- ikT *(l -Af A( F)(r,Y = \n(x /x))) , 



(9) 



where r — \ f\ is the dipole transverse size. In the large- N c limit, one has the following relation between the adjoint 
and fundamental dipoles, 

tf A (r,Y) = Vf F (r,Y)-M${r,Y). (10) 

In the large N c limit, the coupled JIMWLK equations are simplified to the Balitsky-Kovchegov (BK) equation [H|-[l8|, 
a closed-form equation for the rapidity evolution of the dipole amplitude in which both linear radiative processes and 
non-linear recombination effects are systematically incorporated. The running-coupling BK (rcBK) equation has the 
following simple form: 
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where the evolution kernel K run using Balitsky's prescription [l6| for the running coupling is defined as, 
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with r*2 = r — r\ . The only external input for the rcBK non-linear equation is the initial condition for the evolution 
which is taken to have the following form motivated by McLerran-Venugopalan (MV) model @], 
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where the infrared scale is taken A = 0.241 GeV and the onset of small-x evolution is assumed to be at xo = 0.01 [191 ] . 
The only free parameters in the above are 7 and the initial saturation scale Qos (probed by quarks), with a notation 
s = p and, A for a proton and nuclear target, respectively. The initial saturation scale of proton Ql p ~ 0.168 GeV 2 
with the corresponding 7 ~ 1.119 were extracted from a global fit to proton structure functions in DIS in the small-x 
region [li| and single inclusive hadron data in pp collisions at RHIC and the LHC [1, [32|-[34| . Note that the current 
HERA data alone is not enough to uniquely fix the values of Qq p and 7 [l9j]. The recent LHC data, however, seems 
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to indicate that 7 > 1 is preferable [1] . We will consider the uncertainties coming from our freedom to choose among 
different parameter sets for the rcBK description of the proton. 

Notice that in the rcBK equation Eq. (fTTj) . the impact-parameter dependence of the collisions was ignored. Solving 
the rcBK equation in the presence of the impact-parameter is still open problem [35|. However, for the minimum-bias 
analysis considered here this may not be important. Then, the initial saturation scale of a nucleus Qqa should be 
considered as an impact-parameter averaged value and it is extracted from the minimum-bias data. For the minimum- 
bias collisions, one may assume that the initial saturation scale of a nucleus with atomic mass number A, scales linearly 
with A 1 / 3 0], namely we have Qqa = cA 1 ^ 3 Qq p where the parameter c is fixed from a fit to data 1 . In Ref. (34|, it was 
shown that DIS data for heavy nuclear targets can be described with c ~ 0.5. This is consistent with the fact that 
RHIC inclusive hadron production data in minimum-bias deuteron-gold collisions prefers an initial saturation scale 
within Qq A 3 ~ 4 Qq p [H, HH • We will take into account the uncertainties associated to the variation of initial 
saturation scale in the rcBK evolution equation. 



III. OBSERVABLES AND NUMERICAL SETUP 



In this paper, we only consider observables which are defined as a ratio of cross sections. We expect that some of 
the theoretical uncertainties, such as sensitivity to K factors which effectively incorporates the missing higher order 
corrections, will drop out in this way. Therefore, we take K = I throughout this paper. We start by considering the 
ratio of inclusive prompt photon to the neutral pion production, defined as 



^inclusive dN^ ^(Pt) ^ d]ypA—>h(pT)X 
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where the cross-section for the single inclusive prompt photon and hadron production in pA and pp collisions are 
given in Eqs. (I5I7[) . 

In order to investigate the azimuthal angle correlations between the produced prompt photon and hadron, we 
calculate the coincidence probability. In the contrast to a more symmetric production like dihadron, for the photon- 
hadron production we have freedom to select the trigger particle to be a produced prompt photon or a hadron. We 
consider here both cases. In a case that the trigger particle is a prompt photon, the coincidence probability is defined as 
CPh(A(j)) — N^ a " (Acj)) /A^photon, where N^ a "(A(j)) is the yield of photon-hadron pair production including a associated 
hadron (neutral pion) with a transverse momentum p\ s with a trigger (leading) prompt-photon with transverse 
momentum ]P T L and the azimuthal angle between them A(j). In the same fashion, one can define the coincidence 

probability with hadron-triggered as CP 7 (A</>) = A r P air (A</>)/A^h a dron where -/VP air (A^>) is the yield of photon-hadron 
pair including a associated prompt photon and a trigger hadron (neutral pion) with transverse momentum denoted 
by pip s and p\ L , respectively, 



2tt f drPr? f rfJ n /i^^!*i 
CP^{A^;pi. s ,p TL ;T h ,T] h ) = — — —3— , (16) 

J Pt,l d 2 p T h d Vh 

where the integrals are performed within given momenta bins denoted by subscript p T L ,p T s , p T L and p T s . The 
yields in the above expersion are defined in Eqs. (jll4l7p . Similar to the dihadron correlation measurements at RHIC 
[23j, in the definition of the coincidence probability, we follow a convention that a leading or trigger particle has 
transverse momentum larger than an associated particle. Later, we will also study, the implication of different 
kinematics definitions for the trigger particle in 7 — ft, correlations. 

In nuclear collisions, nuclear effects on particle production may be evaluated in terms of ratios of particle yields in 
pA and pp collisions (scaled with a proper normalization) , the so-called nuclear modification factor R p a ■ The nuclear 



Note that a different A-dependence of the nuclear saturation scale with a pre-factor fitted to the HERA data, numerically leads to a 
very similar relation between the proton and nuclear saturation scale [36ll . 
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modification factor for semi-inclusive photon-hadron production is defined as, 

pA t t ^7 ^ N co ii d 2 pr h (Ppr 1 dijy drjh d 2 pr h d 2 pT 7 dr/j drjh 

h h 1 dN pA ^ h( - P 'r')l^T) X dN PP ~^ h( - P T^(PT) X 

R-pA PTiPTiV-yiVh N ^ dptjr dp^, di] 7 dr] h ^ dp\ dp\ dr] 7 dr) h ' ^ ^ 

where the photon-hadron yield in high-energy pA and pp collisions is given in Eq. ([1}. In Eq. (|18[) . the integrals over 
the angles were performed. The normalization constant N co u is the number of binary proton-nucleus collisions. We 
take N co u = 3.6 and 7.4 at yfs = 0.2 and 8.8 TeV, respectively, in pA collisions [37] . Notice that in our approach 
N co u is taken from outset and one should take into account possible discrepancy between our assumed normalization 
N co u and the experimentally measured value for N co u by rescaling our curves. 

We will use the NLO MSTW 2008 PDFs [H and the NLO KKP FFsJp] for neutral pion. For the photon 
fragmentation function, we will use the full leading log parametrization [2^, l40j|. We assume the factorization scale Q 
in the FFs and the PDFs to be equal and its value is taken to be p\ and pj. for inclusive (and semi-inclusive) hadron 
and prompt photon production, respectively. 



IV. MAIN RESULTS AND PREDICTIONS 



In Fig. Q] (right), we show the ratio of single inclusive prompt photon to neutral pion (tt°) production defined via 
Eq. (|14[) for 77^ = fy 7 = 77 and = "p\ = pr at the LHC energy = 8.8 TeV as a function of transverse momentum 
Pt in minimum bias pp and pA collisions at different rapidities rj. It is seen that at the LHC, the ratio j mcluslve [ s 
smaller than one for a large range of rapidities. In Fig. [T] (left), we compare the ratio ry mclusme / ^ a t a fixed rapidity 
r/ = 3 but different energies. In our approach, a fast valence quark from the projectile proton radiates a photon before 
and after multiply interaction on the color-glass-condensate target [l2j . The prompt photon can be mainly produced 
by quark (at the leading log approximation), while pions can be produced by both projectile gluons and quarks, see 
Eqs. (I5I7[) . At the LHC energies at around midrapidity, gluons dominate and consequently the pion production rate 
is higher than prompt photon while for forward collisions and hight pt we have x\ — > 1, therefore projectile quarks 
distributions enhance and consequently the prompt photon production rate grows with increasing rapidity. This can 
be seen from Fig.[T] namely the ratio Y ncluslve /ir° increases with rapidity and transverse momentum while it decreases 
with energy. Note that in our picture, the description of the target appears via the dipole-target forward scattering 
amplitude and it numerically drops out in the ratio, and as a consequence the ratio ^ndustve j n o j g a pp r0 ximately 
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FIG. 2: The photon-hadron (7 — ty°) azimuthal correlation (the coincidence probability) CPh(A<^>) and CP 7 (A<^>) defined in 
Eqs. (| 151 lfcip in minimum- bias (Mini-Bias) pA and pp collisions at RHIC y/S = 0.2 TeV at different rapidities for the produced 
hadron r/h and inclusive prompt photon r] 1 . In the plot, the values of transverse momenta bins of the associated (and leading) 
neutral pion Pt,s (and Pt,l) an d leading (and associated) prompt photon pip L (and pip s ) are given. 



equal for pp and p A collisions at high px and is not sensitive to the initial saturation scale 2 . 

Next, we study the azimuthal angle correlation of photon-hadron (7 — ir°) production by computing the coincidence 
probability defined in Eq. (|15I16I) . In Figs. (|2I3|) we show the coincidence probability for photon-hadron production 
at RHIC and the LHC energy at various kinematics obtained by solutions of the rcBK evolution equation (ITT1) with 
a initial saturation scale for proton Q\ p = 0.168 GeV 2 and for a nucleus within Qqa = 3-^4 Q\ v (corresponding to 
the band) . It is seen that the away-side correlation has a double or single peak structure depending on the definition 
of the trigger (or the leading particle) and kinematics. Namely, if the leading particle is selected a prompt photon 
with pip > pip, then the corresponding coincidence probability CPh(A0) defined via Eq. (|15[) . can have a double peak 
structure at Atfi — it. But if the leading particle is selected to be a hadron with p\ > p T , then the coincidence 



2 Note that in calculation of the ratio of ■y mclualve J ^0 we ignored the inelastic contributions in both inclusive prompt photon and hadron 
production cross-sections assuming that higher order terms will be canceled out in the ratio. 
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FIG. 3: The 7 — it coincidence probability CPh(A(f)) and CPj(Acj>) in minimum-bias pA and pp collisions at the LHC 
y/S = 8.8 TeV at r\h = t? 7 = 3 for two bins of transverse momenta of the associated (and leading) neutral pion s (and l) 
and leading (and associated) prompt photon p^ L (and p^ s ). 



probability CP 7 (Acf>) defined via Eq. p^|) . has a single peak structure at A<fi — ir. In order to understand this 
phenomenon, first note that the cross section of semi inclusive photon-hadron production in quark-nucleus collisions 
given by Eq. ([IJ, becomes zero for: 

PT = \It+Pt~<\ =0. (19) 

This is simply because if the projectile parton is already without any transverse momentum, the production rate 
of photon-hadron should go to zero and off-shell photon remains as part of projectile hadron wavefunction. In 
other words, in order the higher Fock components of projectile hadron wavefunction to be resolved and a photon 
to be radiated, the projectile quark should interact with small-x target via exchanging transverse momentum. The 
necessary kinematics for having a local minimum for the cross-section of photon-hadron production can be readily 
obtained from Eq. (fH))) by using relations given in Eq. ([3]), namely It = p T /Zf and the fact that for the fragmentation 
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FIG. 4: The photon-hadron correlation P(A(f>) defined in Eq. (|22[) in minimum-bias pA collisions at the LHC (right) and 
RHIC (left) at forward rapidity for various transverse momenta of produced prompt photon pj and hadron pj. corresponding 
to two cases of zt > 1 and zt < 1. 



fraction we have z 1 J Lln < Zf < 1. Therefore we obtain, 

z T = V \ < 1, (20) 
Pt 

P V ^ ] < I- (21) 

Note that in our approach, the projectile is treated in the collinear factorization [Hj. Therefore, radiation of photon 
from quark at this level has the standard features of pQCD, including the back-to-back correlation in the transverse 
momentum. Moreover, due to multiple scatterings with target, the cross-section of photon-hadron production should 
have a local minimum for the back-to-back production provided the kinematics conditions given in Eqs. (|20I21|) are 
satisfied. However, because of convolution with fragmentation and parton distribution functions, the local minimum 
will not be zero but gets smeared out. On the other hand, the product oip T Np(pT, x g ) in Eq. (fTJ) has a maximum when 
the transverse momentum px approaches the saturation scale. As a result, a double peak structure appears for the 
away-side correlation. Note that the integrand in the coincidence probability samples smaller transverse momentum 
for photon than hadron for the same reason we already mentioned, namely a photon can be produced if the parton 
already acquired a transverse momentum impulse. In the case that the trigger particle is selected to be a hadron 
rather than a prompt-photon, one should perform the integral over the transverse momentum of the hadron which is 
larger than transverse momentum of photon zt > 1, violating the condition given in Eq. (|20[) . and consequently the 
local minimum at A<j) = tt is washed away and as a result the double peak structure will be fused to a single peak. 
This can be clearly seen in Figs. (12 Kill at both RHIC and the LHC. 

In order to further investigate the consequences of the conditions given in Eqs (|2QI21j) . let us defined the azimuthal 
correlation in the following form [l^ . 

j p A-yfc(j4) 7 (pJ) X j p A->ft(p£) 7 (pJ) X 

P(A^) = — [A<£]/— [A(/> = A<t> c ], (22) 

cPbt p T dp!^ p T dp T dr] 7 drjh d(f> d 2 bt p T dp T p T dp T drj 7 drjh defy 
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FIG. 5: The 7 — n° azimuthal correlation CPh{A(f>) denned in Eq. (|15[) in minimum-bias pA collisions at forward rapidity 
rjh = ^7 = 3 for various energies at RHIC and the LHC. 



which has the meaning of the probability of the semi-inclusive photon-hadron pair production at a certain kinematics 
and angle A<f>, triggering the same production with the same kinematics at a fixed reference angle A<f> c = ir/2. The 
correlation defined in Eq. (|2"21 may be more challenging to measure compared to the coincidence probability defined 
in Eq. (|15I16I) , due to the so-called underlying event dependence, but it its free from the extra integrals over transverse 
momenta and this facilitates to clearly examine the conditions in Eqs. (12012 1[) . In a sense the correlation defined in 
Eq. (|22|) can be considered as a snap shot of the integrand in the coincidence probability defined in Eq. (|15|16|) . 

In Fig. |H we show the photon-hadron correlation P(A(f) defined in Eq. (|22|) at forward rapidity r)h = 77-y — 3 for 
various transverse momenta of produced prompt photon p T , and hadron pfp at RHIC and the LHC for minimum bias 
pA collisions. The initial saturation scale for proton Q^ p = 0.168 GeV 2 and nuclei Qqa — 3Qq p are fixed for all curves. 
It is clearly seen that the photon-hadron away-side correlations can have a double-peak structure both at RHIC and 
the LHC for the kinematics satisfying the conditions in Eqs. (|2QI21j) . and the away-side double-peak correlations will 
evolve to a single peak structure for kinematics outside of region defined by Eqs. ()20l2ip . 

In high-energy collisions, the produced parton on average have intrinsic transverse momentum of order of the 
saturation scale. By increasing the energy or density or decreasing the transverse momentum of the probe, the 
saturation scale Q s increases and consequently this washes away the intrinsic back-to-back correlations and the away- 
side correlation is suppressed. Numerically, a bigger saturation scale, pushes the unintegrated gluon density profile 
to larger transverse momentum. As a result, the single inclusive production (either hadron or prompt photon) cross 
section (the denumenator in the coincidence probability) is enhanced, while the two-particle correlated cross section 
Eq. (dJ is suppressed by a larger saturation scale Q s . Therefore, the coincidence probability defined in Eqs. (|15ll6p . 
decreases with increasing the saturation scale and we expect that the photon-hadron away-side correlation at the LHC 
to be smaller than RHIC (at the same rapidity and transverse momenta of associated and leading particle). This 
can be clearly seen in Fig. [5] where we compare the coincidence probability CPh{A<fi) obtained at various energies. 
Moreover, the saturation scale grows with density, therefore the away-side correlations in pA collisions should be more 
suppressed compared to pp collisions at the same kinematics, see Figs. (|2I3[) . 

It is seen from Figs. (I2I3I4I) that generally at a fixed rapidity and energy, the suppression of away-side 7— h correlation 
is larger for a case that Zt > 1- This effect can be traced back to the fact that 7 — h pairs with zt > 1 probe lower 
a; s -region compared to the cases that Zt < 1. This can be understood by rewriting the definition of x g in Eq. ([3]) 

which appears in the unintegrated gluon density in term of xt, namely x g — -^j= (^e~' n ~< + |y e^ 71 ^ . Therefore, 7-/1 
pairs production with zt > 1 have a lower x g and consequently the suppression due to saturation will be larger. 

In Fig. H we show CP h (A(f>) and CF 7 (A0) at the LHC s/S = 8.8 TeV at forward rapidity r\ h = ^ = 3 in 
minimum-bias pA and pp collisions for two bins of transverse momenta of associated prompt photon p T s and hadron 
Pj. s , and the corresponding leading hadron pip L and prompt photon p^ L . Namely in top and lower panel we 
performed the integral for the associated particle within [1,2] GeV and [1,6] GeV (and for the corresponding leading 
particle within [2, 20] GeV and [6, 20] GeV), respectively. The correlation signal enhances by increasing the transverse 
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FIG. 6: The 7 — 7r° azimuthal correlation CPh(Acj)) defined in Eq. (|15p in minimum-bias pA and pp collisions at the LHC 
y/S = 8.8 TeV at different rapidities of the produced hadron rjh and prompt photon rj^. 



momenta interval of associated particle. This is simply because in Eqs. (|15I16[) . by construction, the integrals over 
the leading particle is mainly canceled out in the ratio, and the correlation becomes proportional to the integral over 
the associated particle. For higher transverse momenta bins, the saturation scale is smaller and the back-to-back 
correlation is restored. Notice that since the rcBK evolution solution is not reliable at high transverse momentum 
we had to impose upper limit cut for the integrals over transverse momenta of the leading particle in Eq. (|15I16I) . 
However, the cross-sections drop so fast with transverse momentum at forward rapidities that this upper cutoff should 
not make a big difference. 

The azimuthal correlation CPh(A<fi) defined in Eq. (TT5l) is generally bigger than the corresponding correlation 
CP 7 (A0) defined in Eq. (fTB|) . at the same kinematics. This is because when the trigger particle is taken a prompt 
photon, in Eq. (|T5|). the electromagnetic coupling a em drops out in the ratio of two cross sections and that enhances 
the signal compared to the case that the trigger particle is selected to be a hadron. This can be seen in Figs. (|2I3I) . 

As we already pointed out, the double peak structure for the photon-hadron coincidence probability CPh{A-4>) at 
A(f> « 7T is due to the interplay between a local minimum for the cross-section at px ~ and two maxima for the 
cross-section when px ~ Q s - The double-peak structure can be stretched out and becomes more pronounced by 
measuring the associated hadron at about or higher rapidity than the trigger prompt photon, i.e. r\h > V-y- This is 
due to the fact that because of kinematic limit for more forward production, the integrand of the associated hadron 
in CPh ( A0) is relatively shifted to lower transverse momentum and consequently the conditions for local minimum in 
Eqs. (|20I21[) are satisfied while at the same time, the saturation scale increases for more forward production leading to 
an enhancement of the two local maxima. In Figs. (|2l6p . we show this effect by comparing the azimuthal correlations 
at different rapidities rjh and ry 7 at RHIC and the LHC. 

Although the main features of the photon-hadron correlations, e.g. the double or single peak structure and decor- 
relation with energy /rapidity, density and transverse momentum seem to be robust and understandable due to the 
non-linear gluon saturation dynamics, there is some uncertainties on the magnitude of the correlation obtained in our 
approach. These uncertainties are due to the fact that with available worldwide small-x experimental data it is not 
yet possible to uni quely f ix the parameters of the rcBK evolution equation and the initial saturation scale of proton 
and nucleus [8l. [lil, l32T[34j . To highlight our main uncertainties, in Fig.[7]we show CP/j(A</>) for minimum-bias pA and 
pp collisions at forward rapidity at RHIC and the LHC, with two different initial saturation scale of proton, namely 
Qq p = 0.168 and 0.2 GeV 2 corresponding to 7 = 1.19 and 7 = 1 in Eq. (fHij) respectively [l9j], and the initial saturation 
scale of nuclei (gold and lead) within Q\ A = 3 ~ AQq . Although both values of Qo p (or 7) are extracted from a 
fit to HERA data on the proton target at small-x [19| (with a good x 2 )> the recent LHC data seems to favor the 
parameter set with 7 > 1 (or lower value for Qo p ) Q. The uncertainties in the initial scale for proton will bring even 
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A<|> A$ 

FIG. 7: The 7 — 7r° azimuthal correlation CPh(A0) in minimum- bias pA and pp collisions at forward rapidity at RHIC and 
the LHC. The curves are obtained by the rcBK equation with two different initial saturation scale of proton Qq p = 0.168 and 
0.2 GeV 2 and the corresponding initial saturation scale of the nucleus within Qf )A = 3 + 4Qq p . 



larger uncertainties in determining the parameters of the rcBK equation for the case of nuclear target 3 . Therefore, the 
upcoming LHC data on pA collisions can provide crucial complementary constrain on the rcBK evolution equation 
and small-x physics. For other measurements sensitive to the saturation physics at the LHC, see Refs. foL [T21 132^. 

Note that the semi-inclusive photon-hadron cross-section in Eq. (JTJ has collinear singularity. Therefore, one should 
first treat the collinear singularity for the near-side jet Aip ss in a same fashion as was done for the inclusive prompt 
photon production in Eq. Q by introducing the quark-photon fragmentation function. Therefore, our results at 
near-side Acj) « should be less reliable. However, one should bear in mind that the integrand in the azimuthal 
correlation generally samples lower transverse momenta for the away-side correlations A(j) rs tt than for near-side ones 
A<p ~ 0. Therefore, here we only focused on the away-side correlations which is a sensitive probe of small-x physics 
and gluon saturation. 

In Fig. [8l we show the nuclear modification factor R p ^ for semi-inclusive photon-hadron pair production defined in 

Eq. (|17p as a function of Aip at the LHC energy y/S = 8.8 TeV and rjh = ry 7 = 3 for two different transverse momenta 
bins of produced prompt photon p7p and hadron p^ (the integral is performed over the given interval of transverse 
momenta). Similar to previous plots, the band (CGC-rcBK-av) in Fig. [8] comes from the rcBK solutions incorporating 
the uncertainties associated to a variation of the initial saturation scale of the nucleus in a range consistent with 
previous studies of DIS structure functions as well as particle production in minimum-bias pp, pA and AA collisions 
in the CGC formalism. One may therefore expect that the possible effects of fluctuations on particle production is 
effectively contained in our error band. The away-side nuclear modification R p J^ at Atfi ss it is dramatically suppressed 
with a lower peak structure when the transverse momentum bin of the produced prompt photon is larger than hadron 
Zt < 1- This is fully in accordance with the photon-hadron decorrelation in pA compared to pp collisions, and 
conditions given in Eqs. (|20I21[) for the existence of the local minimum for the away-side photon-hadron production. 
Note that the sensitivity to the transverse momenta or the ratio zt only manifests itself at around Acj> ~ 7r. 

Finally, in Fig. [§1 we show the two-dimensional nuclear modification factor for semi-inclusive photon-hadron 
production defined in Eq. (|18[) as a function of transverse momentum of produced prompt photon pj, and hadron 
p\ at RHIC \[S = 0.2 TeV at % = T} 7 = 4 (top panel) and at the LHC \/S = 8.8 TeV at r) h = ^ = 3 (lower 
panel). The area between two surfaces in Fig. [9j similar to Fig. [8] (the band labeled by CGC-rcBK-av) shows the 
uncertainties associated to the variation of the initial saturation scale of the nucleus. It is seen that at the LHC 
energy y/S = 8.8 TeV, the nuclear modification factor i?^ is more suppressed compared to RHIC and also is more 
flat. We recall that the semi-inclusive photon-hadron cross-section Eq. |T]) is not equal to the product of cross-sections 



3 This is partly due to the fact that solution of the rcBK equation in the presence of impact-parameter is not yet available. 
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FIG. 



The nuclear modification factor for semi-inclusive photon- hadron (7 — ty°) production defined in Eq. (|17[) as a 



function of A(f> at the LHC in minimum-bias pA collisions at forward rapidity r\n = V-y = 3 for two different bins of transverse 
momenta of produced prompt photon an d hadron pj., namely zt < 1 (right) and zt > 1 (left). The band (CGC-rcBK-av) 
incorporates the uncertainties due to variation of the initial saturation scale in the rcBK evolution equation. 



of single inclusive prompt photon and hadron production given in Eq. (0} and Eq. J7J). Note that in Ref. [lj] it was 
shown that at RHIC for the single inclusive prompt photon production, a good portion of the suppression at forward 
rapidities is due to the projectile being a deuteron rather than a proton. Here for a comparison with pA run at the 
LHC energy and in order to discard possible suppression associated to isospin effect [l2j . we have only considered 
proton- nucleus collisions at the RHIC energy which can be also useful for the future pA run at RHIC. We check that 
similar to Ref. [12| , discarding the fragmentation photon contribution from the cross-section, will not affect our results 
for R!l\ significantly. Nevertheless, a detailed study of the semi-inclusive photon-hadron production in the presence 
of isolation cut is beyond the scope of the current paper. 



V. SUMMARY 



We have investigated semi-inclusive prompt photon-hadron production in high-energy pp and pA collisions within 
the CGC framework by using the running-coupling BK equation. We provided detailed predictions for the coincidence 
probability of photon-hadron correlations and showed that such correlations exhibit novel feature, namely the away- 
side correlations can have a double or single peak structure depending on the trigger particle selection and kinematics. 
The correlations have a double-peak structure by selecting 7 — /i pairs within the kinematics region satisfying the 
conditions in Eqs. (|20I21I) . and the away-side double-peak correlations will evolve to a single peak structure for 
kinematics outside of that region. We showed that this feature can be understood by QCD saturation dynamics. 
The double-peak structure for the azimuthal correlations has been also recently reported for other electromagnetic 
probe, namely the Drell-Yan Lepton-Pair- Jet correlation in pA collisions [4l| while it is absent for dihadron production 
[131 [23j . The decorrelation of the away-side photon-hadron production with energy, rapidity, density and transverse 
momentum of the probe is very similar to the dihadron production in pA collisions and can be understood in the 
CGC framework. If experimentally confirmed, this will provide a significant evidence in favor of the universality of 
particle production in the QCD saturation picture at high-energy. 

In a sense, the double-peak structure for 7 — h correlations resembles the long-range azimuthal correlations for the 
produced charged hadron pairs, observed in high-multiplicity events in pp collisions at the LHC, the so-called ridge 
phenomenon [42[. Although, the ridge is a feature on a near-side Acj> ~ of the two particle correlations, while the 
7 — h double-peak structure is a away-side feature. In both cases, a second local maximum occurs because of angular 
collimation due to the presence of the saturation scale in the system, and the effect shows up within a kinematics 
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FIG. 9: The nuclear modification factor R p ^ for semi-inclusive 7 — n° production defined in Eq. (|18[1 as a function of transverse 
momentum of produced prompt photon p^, and hadron p\ in minimum-bias pA collisions at RHIC yfS = 0.2 TeV at r) h — -q-y = 4 
(top panel) and the LHC y/S = 8.8 TeV at r)h = Vj = 3 (lower panel). Two surfaces are obtained from the solutions of the 
rcBK evolution equation with two different initial saturation scale of the nucleus, see the text for the details. 



window which is dictated by the saturation scale 4 [l(J HH, IHI • Similar to the ridge, the double-peak structure here 
can survive up to rather large rapidity (see Fig. [5]), and in both cases, one expects that the same mechanism to be 
responsible for the self-deconstruction of the effect namely decorrelation at very high-energy [43j . 

We also showed that the ratio zt = Px/Pt 1S a sensitive parameter to the saturation region and controls the 
away-side 7 — h suppression in high-energy pp and pA collisions. 

We studied the ratio of single inclusive prompt photon to hadron production ^ mcluslvc / 7r ° m pp anc } pA collisions at 
RHIC and the LHC at various rapidities. We found that the ratio ^ mclusive / 7r ° i s very similar for high-energy pp and 
pA collisions at forward rapidities at high transverse momentum, and it increases with rapidity while it decreases with 
energy. We also provided predictions for the nuclear modification factor for the semi-inclusive photon-hadron pair 
production R p \ in pA collisions at RHIC and the LHC at forward rapidities. We showed that the two-dimensional 

RpA 1S generally more flat at the LHC compared to RHIC at forward rapidities. We found that the suppression of 
the nuclear modification factor for semi-inclusive photon-hadron production is comparable to that for single inclusive 
hadron [32j and prompt photons [l2| production in pA collisions at forward rapidities. 



4 The author thanks Raju Venugopalan for pointing out the possible similarity between these two phenomena. 
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